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ORBITAL PARAMETERS 
FOR THE OGO-E 

GEO-CORONAL HYDROGEN EXPERIMENT 
ABSTRACT 

This document presents and applies methods for obtaining 
certain necessary orbital and spacecraft parameters for meaning- 
ful data analysis of the geo-coronal hydrogen measurement 
experiment (Experiment No. 22) on the fifth in Orbiting Geophysical 
Observatory Series, namely the OGO-E. 

Techniques are presented for determining as a function of 
time: l) orbit inclination with respect to the ecliptic plane, 

2) the angle between the orbit's angular momentum vector and the 
sun's position vector, 3) the half angular dimension of the earth 
viewed from the satellite, 4) the solar array angle and, 5) the 
angle between the projection of the sun on the x-y plane of the 
experiment's coordinate system and the x-axis of that system. 

The calculations were made using a numerical integration 
general purpose interplanetary trajectory program (ITEM). Conse- 
quently, the changes in the orbit due to perturbations are reflected 
in the calculations of the quantities mentioned above. 

The results of computer runs are presented in this document. 
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NOMENCLATURE 


a 


Semi -major axis. 


a e 


Radius of the earth. 


Ce 


e 


The angle between the OEEP xg-axis and the 
projection of the sun vector on the OPEP Xg-yg 
plane. (See Figure 2.) 

Eccentricity. 


i 


Inclination. 


i, 3, k 



:L b , 'V \ 

3 E * kg 



M 

N 

r 

S’ 




Unit vectors along the x, y, z axes of the equa- 
torial inertial coordinate system. 

Unit vectors along the x 1 , y', z* axes of the 
inertial ecliptic coordinate system. 

Unit vectors along the xg, y-^, z^ axes of the 
satellite fixed coordinate system. 

Unit vectors along the xg, yg, zg axes of the 
OPEP coordinate system. 

Unit vectors along the Xp, yp, Zp axes of the 
solar array coordinate system. 

Mean anomaly. 

Unit angular momentum vector. 

Position vector of the satellite. 

The angle between normal to the orbit and the 
position vector of the sun relative to the earth. 

Sun vector projection on the OPEP Xg- y E 
plane. 

Unit vector from earth to sun in the x, y, z 
coordinate system. 

Axes of the inertial equatorial coordinate system. 
Axes of the inertial ecliptic coordinate system. 


x, y, z 
x*, y*, z' 


Vll 



x b> y b > z b 

X E J y E > Z E 

x , y , z 
P P P 

Greek Symbols 
6 


? 

^E 

0) 

Q 


Axes of the satellite fixed coordinate system. 
Axes of the OPEP coordinate system. 

Axes of the solar array system. 


Half angular dimension of the earth viewed from 
the satellite. 

The mean obliquity = 23?V+37 for 19^6 (Ref. 2 ). 

The orbit inclination with respect to the ecliptic 
plane . 

The solar array angle (Fig. 2). 

The OPEP angle, the angle between the Xg and x^ 
axes . 

Argument of perigee. 

Right ascension of the ascending node. 


Matrices 

A 

B 

C 

X 


X 


b 


X E 


A transformation matrix defined by equation 8. 

A transformation matrix defined by equation 9 - 
= B Aj see equation 10. 

A vector with components in the inertial equatorial 
coordinate system. 

A vector with components in the satellite fixed 
coordinate system. 

A vector with components in the OPEP coordinate 
system. 



Position vector of the sun in the OPEP system. 


Position vector of the sun in the inertial 
equatorial coordinate system. 
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ORBITAL PARAMETERS 
FOR THE OGO-E 

GEO-CORONAL HYDROGEN EXPERIMENT 
INTRODUCTION 

The Orbiting Geophysical Observatory (OGO) satellites 
form a family of satellites each designed to carry up to 
fifty experiments. The OGO-E is one in this series. The 
discussion here concerns an experiment (No. 22) on the OGO-E 
which will make optical measurements of the thickness of the 
geo-coronal hydrogen. It is theorized that hydrogen is 
escaping from the earth. The purpose of the experiment is 
to determine the gradient of the hydrogen population around 
the earth and also to investigate the interaction caused by 
extra-terrestrial effects . 

This document gives methods for determining the following 
quantities useful in post flight data analysis: 

1. The angle £ between the ecliptic plane and the 
orbit plane (Fig. l) . 

2 . The angle S * between the normal to the orbit plane 
and a vector from the earth to the sun (Fig. l) . 
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3. The half angular dimension 6 of the earth as viewed 
from the satellite as a function of time . 

4. The solar array angle cp between the normal to the 
sun-lit side of solar array and the y^-axis of the 
satellite fixed coordinate system (Fig. 2). 
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5. The angle C E between the projection of the sun 
vector on the x^,- y^ plane of the orbit plane 
experiment package (OPEP) coordinate system and 
the Xg-axis of that coordinate system. The OPEP 
coordinate system, x^,, y E , z E is shown on Figure 2. 



Figure 2. 

Angles and C £ 
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I. ASSUMPTIONS 


The OGO-E will be placed into a highly eccentric orbit 
with an apogee height approximately 148,000.0 km. The injection 
conditions for OGO-E are assumed to be the same as the nominal 
conditions for 0G0-A which was launched September 4, 19 64 (Ref. l) . 


Table I. 

OGO-E Injection Conditions 


Geocentric latitude 

= 

209744573 s 

Longitude 

= 

111911923 E 

Height 

= 

279.2517 km 

Speed 

= 

IO.716286 km per sec 

Azimuth 

= 

669445986 from north 

Flight path angle 

= 

19452111 


The osculating classical orbital elements associated with 
these injection conditions are shown in Table II. 
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Table II 


Orbital Elements 


Semi -major axis, a = 12. 51*1-687 earth radii 

Eccentricity, e = .91666199 

Inclination, i = 30?910^96 

Argument of perigee, tw = 313?657536 

Right ascension of the ascending node, 0 = 

Mean anomaly, M = 8?7l60632 

Epoch = August 15, 1966 at 5^ 30 m U.T. 


195?7*K)31 
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II . MATHEMATICAL MODEL 


A. Orbit Plane Geometry 

The angles § and S 1 are related to the orbit in space 
and the sun. These angles are considered as orbit related 
because the position of the spacecraft in the orbit does not 
enter into the calculations . 

1. The Inclination with respect to the Ecliptic Plane I. 

The first problem is to determine the angle between the 
normal to the ecliptic plane and the normal to the orbit plane, 
i.e. the inclination of the orbit with respect to the ecliptic 
plane. The angle can simply be written as: 


§ = cos 1 (N • k ’ ) • 


( 1 ) 


where N = Unit angular momentum vector which is 
normal to the orbit. 

k' = Unit vector normal to the ecliptic plane, 


k' 


= - j sin e + k cos e 


e = 23?44371, mean ecliptic angle for 19 66 (Ref. 2) 

j, k = unit vectors along the y and z axes of the 
equatorial inertial coordinate system. 
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2. Angular Momentum - Sun Angle, S'. 

The angle S ’ is defined by: 

S' = cos -1 (I • U gE ) (2) 

where Ug E = The unit position vector of the sum measured 
from earth. The coordinates of the sun are 
taken from an ephemeris. 

The use of N in equations (l) and (2) takes into account 
the motion of the orbit due to perturbations when used in a 
numerical integration scheme where for each integration step 
N is computed. Consequently, the motion N reflects both the 
secular and periodic perturbations. For a single orbit there 
is little change in N; hence, it can be stated that the angles I 
and S’ are dependent on the motion of the plane of the orbit. 

The angle S' is also dependent on the position of the sun. 

B. Spacecraft Angles 

1. Half Angular Dimension of the Earth, 6. 

The angle, 6 is easily determined from the following 
expression: 

a 

6 = sin" 1 0° ^ 6 £ 90 ° ( 3 ) 

M 

where a g = 6378 . 16 5 km radius of the earth (Ref. 3 ). 
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2. The Solar Array Angle, 9 . 

P 

The solar array angle, cp^, is based solely on the position 
of the solar array with respect to the main structure on the 
satellite. The angle 9 is measured counter-clockwise from 
the y^-axis to the y^-axis, as observed from the positive 
x^-axis (Fig. 2). The planned range for the variation of 
9 for OGO-E is ^ ^ 9 s ^II, (Refs. 1*-, 5). Therefore, 9 may • 
be found from the formula 


P 


cos -1 (jp 




no 


where 


Unit vector normal to solar array surface 
(sunlit side). 


*b 



r is the satellite position vector. 


3. OBEP Xg-Axis - Sun Projection Angle, C E . 

The angle C_ is somewhat more complicated in that two 

Ej 

coordinate transformations are required involving three 
coordinate systems. The three coordinate systems in question 
are: 


a). X 


— - 

x 

y 

z 


inertial equatorial system, 
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where x, y, z form the right-handed coordinate system in 
which x points to the vernal point, z points to the north 
pole and y completes the set. 

, the satellite fixed coordinate system. 



where x^, y^, z^ form a right-handed coordinate system 
with the origin at the center of gravity of the satellite 
in which z^ points to the center of the earth from the 
satellite, x^ lies parallel to the solar array axis and 
y^ is orthogonal to x^ and z^ and points towards the 
orbit plane experiment package (OPEP) end of the satellite 
(Fig. 2). 


c). Xj 



the OPEP coordinate system. 


where Xg, y E , Zg form a right-handed coordinate system in 
which Zj, is parallel to z^ and the Xg - y £ plane is 
parallel to the x^ - y^ plane and has the angle Yg between 
the x^ and Xg axes as shown in the following figure. 
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In the three coordinate systems described above, let: 


i, 3 , k = unit vectors along x, y, z. 

V V \ = unit vectors along x b , y b , z b . 
i E , I E , kg = unit vectors along Xg, y E , Zg. 
Now, the matrix A is defined by 



V T 

V 1 

V s 

A = 

V 1 

V 7 

V* 


v r 


V E 

matrix B is: 




cos Y e 

sin 


B = 

- B in Y e 

cos 



0 

0 



0 

0 • 
1 


( 8 ) 
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Therefore, 


C 


where 


Cll 

C12 

C 13 

C21 

C22 

C 2 3 

C31 

C32 

C33 



— 


cii = (i b -i) cos Y e + sin Yg 

Cx 2 = (i b *tD cos Yg + sin Yg 

cia = (i b *k) cos Y E + (J b *k) sin Yg 
Cgi = — ( i b * i ) sin Yg + (j b *i) cos Yg 

c 22 = -(i b *3) sin Yg + ( 0 b * J) cos Yg 

c 23 = -(I b -k) sin Yg + ( j b -k) cos Yg 

c 3 i = 

C32 = 

C33 - (k^'k)’ 
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Since the position of the sun must be defined in the OPEP 
coordinate system, equation (7) is used. 



and Xg, yg, Zg are the coordinates of the sun. Hence the 
coordinates of the sun in the OPEP system are: 



In terms of vector notation the sun vector in the OPEP 
system is: 


“ ^E ^g + ^E y Eg + k E Z Eg 


( 12 ) 


13 



The projection of the sun onto the OPEP Xg - y E plane is 
simply: 


3 E^ " i E X E r , + J E y ] 


E, 


( 13 ) 


The angle Cg between the Xg axis and the projection of the sun 
on the OPEP Xg - y £ plane defined by: 

S T 


C E = cos 


E 


Is. 


( 1 ^) 


E 


C. The Computer Program 

The equations in Section II were used with an existing 
general purpose interplanetary trajectory Encke method program 
(ITEM) for use on the IBM 709O and IBM 709*4- computers (Ref. 6). 

The perturbations accounted for in this program include the 
gravitational attraction of the earth, moon, sun, Jupiter, Venus, and 
Mars considered as point masses. Also included are the effects 
of the 2nd, 3rd and 4th zonal harmonics and 2nd tesseral harmonic 
of the earth’s gravitational field. Since perigee height rises 
rapidly, air drag was not considered. The perigee height rise 
for OGO-A is about fourty-one kilometers per orbit. 
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III. RESULTS AND DISCUSSION 


Computer runs to determine §, S', 
functions of time were made using the 
and the following injection times: 


6, and C E as 
formulas of Section II 


August 7, 1966 
August 7, 1966 
August 22 , 1966 
August 22 , 1966 


1H5” 

universal time 

3 h 45” 

universal 

time 


universal 

time 

fl 5 " 

universal 

time 


A. Orbit Plane Geometry 

1 . Inclination with respect to the Ecliptic Plane, §. 

Figures k and 5 show the trend of I for the above injection 
times. This angle exhibits the approximate fourteen day 
lunar oscillation and the bi-annual oscillation due to the 
solar perturbation. The mean increasing trend is due 
primarily to the combined long term effects of the sun and 
moon. 
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ANGLE $ (Degrees) Angle J (Degrees) 



2. Angular Momentum - Sun Angle, S* 

Since S’ is mainly a function of the sun's position, 
the time history of S' shows a yearly periodicity. The short 
term periodic (about fourteen days) lunar perturbation also 
affects S' but these variations are in the order of tenths 
of a degree (about two orders of magnitude less than the 
average value of S') and, consequently, are not perceptible 
on Figures 6 and 7 . 
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ANGLES' (DEGREES) ANGLES' (Degrees) 





B. Spacecraft Angles 

1. Half Angular Dimension of the Earth, 6. 

The angle 6 is nearly periodic since it is a function 
solely of the distance of the spacecraft from the earth. 

The time history for the first 150 hours after injection is 
essenti ally the same for all the injection dates and times 
shown on page 15. See Figure 8. 

After almost one year of flight the time history of 
6 is again shown for 150 hours (Figure 9). These results 
are based on an injection into orbit time of August 7, I966, 
at J* 1 I5 111 universal time. The results shown on Figure 9 
are typical of results for different injection dates and 
times . 
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ANGLE & (DEGREES) ANGLE ^ (DEGREES) 
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FIGURE9 - & Versus Time after Injection 






2. The Solar Array Angle, cp^. 

Figures 10 and 11 show the behavior of the solar array 
angle cp^ for the first 150 hours in orbit for the injection 
into orbit times shown on page 15 • The different time 
histories are due to the different spacecraft-sun geometry. 

Figures 12 and 13 show the behavior of the solar array 
angle cp^ after alm ost a year in orbit. 
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SOLAR ARRAY ANGLE %> ( DEGREES ) 



TIME AFTER INJECTION ( HOURS ) LAUNCH DATE (8-7-66) 


FIGURE 10 - ff> Versus Time after Injection 



FIGURE II- Versus Time after Injection 
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SOLAR ARRAY ANGLE. f f ( DEGREES ) SOLAR ARRAY ANGLE, ft (DEGREES) 




J. OPEP x^-Axis - Sun Projection Angle, Cg. 

Figures Ik and 15 show the behavior of C E for the 
first 150 hours in orbit for the injection into orbit times 
shovn on page 15 . As in the case of the solar array angle 
the difference in the time histories of C £ for different 
injection times is due to the different spacecraft-sun 
geometry. 

Figures 1 6 and 17 show the behavior of after almost 
a year in orbit. 
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ANGLE C.- (DECREES) 
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